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Executive  Summary 


In  the  quest  for  development  of  structural  materials  with  high  temperature  capability  well 
above  1000°C  in  primary  structural  applications  many  systems  have  been  under  considera¬ 
tion.  Of  these,  among  ceramic  systems,  the  directionally  solidified  (DS)  eutectic  ceramics  have 
demonstrated  interesting  potential  for  use  in  the  1200-1500°  range  as  competitors  to  sapphire 
fiber.  A  system  among  these  eutectics  that  has  such  potential  but  also  is  relatively  simple  in 
composition  and  morphology  is  the  DS  eutectic  of  Al203/c-Zr02(Y203). 

In  the  as-solidified  form  this  eutectic  combines  a  topologically  continuous  majority  phase  of 
AI2O3  at  a  volume  fraction  of  c.a.  0.67  with  a  nearly  perfect  texture  of  [0001]  parallel  to  the 
growth  axis  that  encapsulates  a  minority  phase  of  cubic  Zr02  in  a  submicron  fiber  or  platelet 
morphology.  The  low  temperature  fracture  toughness  of  this  eutectic,  at  a  level  of  c.a.  5MPa 
m1/2  is  considered  adequate.  The  behavior  of  primary  importance  of  this  eutectic  for  structural 
applications  is  its  creep  resistance  in  the  1200-1500°  range  within  which  the  morphological 
stability  against  coarsing  is  acceptable.  It  is  the  creep  resistance  of  this  eutecitc  that  has  been 
under  study  in  this  AFOSR-supported  research  program. 

Based  on  detailed  morphological  characterization  of  the  eutectic  and  on  some  associated 
computer  simulations  of  dislocation  cores  as  well,  as  FEM  studies  of  the  internal  stress  distribu¬ 
tions  among  phase  components  under  stress,  the  creep  experiments  in  the  range  of  1200-1520°C 
at  stress  levels  up  to  300MPa  resulted  in  a  creep  model  that  is  in  very  good  agreement  with 
experimental  findings.  The  creep  model  is  based  on  the  response  of  the  topologically  continu¬ 
ous  AI2O3  component.  It  proposes  that  creep  in  the  AI2O3  component  results  entirely  from  the 
diffusion  controlled  climb  of  pyramidal  edge  dislocations  which  are  forced  to  thread  through 
the  sub-micron  c-Zr02  domains  acting  as  dispersoids.  The  theory  which  agrees  well  with  the 
experimental  measurements  of  steady  state  creep  in  the  1200-1520°C  range  predicts  that  in 
this  range  the  eutectics  have  better  creep  resistance  than  pure  sapphire  of  ”c”  axis  orientation, 
provided  that  the  morphology  does  not  coarse.  This  latter,  possibly  adverse,  behavior  needs 
further  study. 

The  exploitation  of  the  very  good  creep  resistance  of  these  eutectics  in  directionally  solidified 
fiber  form  will  need  substantial  additional  engineering  development  before  material  components 
can  be  made  ready  for  structural  application.  Much  earlier  developments  on  monazite  coatings 
to  prevent  surface  degradation  of  eutectic  fibers  in  compositions  should  make  such  engineering 
materials  systems  developments  more  readily  possible. 

The  research  at  the  Massachusetts  Institute  of  Technology  was  aided  importantly  by  the 
collaboration  of  A.  Sayir  of  the  NASA  John  Glenn  Research  Center  who  supplied  all  the  DS  eu¬ 
tectics  and  also  contributed  to  the  interpretation  of  the  creep  results.  Other  financial  assistance 


1  Introduction 


To  reach  beyond  the  capabilities  of  the  present  set  of  high  temperature  superalloys  and  some 
intermetallic  compounds  that  have  service  limitations  at  around  1000°C,  it  is  necessary  to 
consider  other  compounds  such  as  oxides,  carbides,  borides,  etc.  to  attain  service  temperatures 
in  the  1400-1600°C  rangefl].  While  these  are  all  intrinsically  brittle  materials,  having  low 
temperature  brittleness  problems,  their  high  temperature  performance  is  largely  governed  by 
their  creep  resistance  and  fracture  resistance.  It  is  creep  resistance  that  was  of  principal  concern 
in  the  completed  research  program. 

For  some  time,  single  crystal  sapphire  fiber  with  [0001]  axis  orientation  has  been  considered 
as  an  ideal  material  for  high  temperature  application  [2].  In  that  orientation  with  its  principal 
basal  and  prismatic  systems  unstressed,  sapphire  single  crystal  fibers  have  remarkable  creep 
resistance  in  the  1400-1600  °C  temperature  range  where  only  the  pyramidal  slip  system  is 
stressed.  Early  experiments  of  Firestone  and  Heuer  [3]  on  [0001]  axis-oriented-sapphire  in  the 
1600-1800°C  range  produced  evidence  that  such  sapphire  crystals  creep  entirely  by  the  climb  of 
the  l/3<  1101  >  pyramidal  edge  dislocations,  with  no  slip  line  or  stereo-TEM  evidence  of  glide 
of  such  dislocations  on  any  of  the  possible  pyramidal  planes  available.  That  such  dislocations 
are  entirely  sessile  in  glide  has  now  been  established  in  very  recent  MD  simulations  of  the  core 
structures  of  these  dislocations  in  the  present  program,  and  is  presented  in  Section  4.2  [4,  5]. 
The  stress  dependence  of  the  creep  rates  of  sapphire  crystals  of  [0001]  orientation  and  their 
governing  activation  energy  being  of  oxygen  ion  diffusion  in  AI2O3  have  all  been  consistent 
with  a  pure  climb  mode  of  creep  [3].  This  is  extremely  rare  in  the  creep  of  metals,  with  the 
only  similar  occurrence  having  been  reported  by  Edelin  and  Pourier  [6]  in  similarly  oriented 
Mg  crystals. 

While  oriented  sapphire  single  crystal  fibers  have  advantageous  characteristics  and  have 
been  considered  with  diffusion  barrier  coatings  as  reinforcing  elements  in  polycrystalline  AI2O3 
matrixes  for  composite  applications  [7,  8],  the  practice  has  been  costly  and  lacked  flexibility. 
As  an  alternative  to  sapphire-fiber-reinforced  ceramic  composites,  a  series  of  directionally  so¬ 
lidified  advanced  ceramic  eutectics,  consisting  largely  of  an  AI2O3  component  together  with 
compatible  stable  oxides  of  ZrC>2  or  YAG  combine  many  of  the  advantages  of  sapphire  fiber 
with  morphological  stability  at  elevated  temperatures  in  the  1500°C  range  and  relative  ease  of 
production.  Of  these  the  AI2O3,  Zr02  system  eutectic,  with  additional  4.2  mol%  Y2O3  modifi¬ 
cation  to  create  cubic  Zr02,  has  demonstrated  many  attractive  characteristics  [9].  This  system 
was  chosen  for  detailed  study  in  the  present  research  program. 
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2  Experimental  Details 

2.1  Material  and  characterization 

The  eutectics  of  A^Oa/c-ZrC^^Oa)1  have  all  been  produced  by  the  Laser  Heated  Float  Zone 
(LHFZ)  method  at  the  NASA  John  Glenn  Research  Center  in  the  facilities  of  our  principal 
collaborator  A.  Sayir  [9].  Source  rods  were  prepared  using  99.999%  pure  polycrystalline  AI2O3 
powder  (CERAC/pure,  Ceralox  Corp.,  Tucson,  AZ  08676),  and  99.999%  pure  ZrC>2  and  Y2O3 
powders  (Alfa  Aesar).  These  powders  were  blended  in  acetone  for  70  hours.  The  slurry  was 
dried,  mixed  again  and  formed  into  cylindrical  rods  using  cold  isostatic  pressing.  These  rods 
were  then  sintered  at  1500°C  in  air  for  4  hours.  A  LHFZ  process  was  used  to  produce  the 
directional  solidification  of  rods  as  described  previously  [2].  The  rods  were  re-used  as  source 
rods  to  produce  final  samples.  The  objective  of  the  second  solidification  was  to  minimize  the 
entrapment  of  gases  at  the  liquid  solid  interface.  The  directional  solidification  (DS)  rate  was  at 
2mm/min  for  both  processes.  The  final  DS  rods  had  cross-sectional  dimensions  roughly  in  the 
range  of  1-1.5  mm  diameter  changing  from  rod  to  rod  but  showed  only  small  thickness  variations 
along  their  20-22  cm  length  in  any  one  rod.  They  were  composed  of  67%  AI2O3  and  33%  c- 
ZrC>2  by  volume.  Random  sectioning  of  the  rods  always  showed  a  certain  residual  component  of 
Y3AI5O12  (YAG)  at  levels  that  were  considered  to  be  insignificant.  More  significant  however  was 
an  unavoidable  level  of  porosity  along  the  center  line  of  many  samples,  often  with  substantial 
pore  dimensions  to  constitute  super-critical  flaws  for  fracture.  These  pores  always  had  smooth 
surfaces  indicating  they  resulted  from  inadequate  wetting  during  melting  and  solidification  of 
the  initial  charge.  In  a  few  instances  they  resulted  in  fracture  during  loading  of  the  creep 
experiments  where  they  were  given  no  further  attention  other  than  characterization  of  the 
pore  sizes  and  shapes  as  a  feed-back  to  improve  the  production  process.  TEM  specimens 
made  from  randomly  spaced  axial  and  transverse  sections  were  used  for  electron  diffraction 
to  check  orientation  of  components.  These,  as  well  as  a  limited  number  of  X-ray  pole  figure 
determinations  obtained  from  surfaces  of  axial  and  transverse  section  of  the  bars  demonstrated 
that  the  AI2O3  topologically  continuous  phase  had  a  growth  texture  of  [0001]  within  2-3°  parallel 
to  the  specimen  axis  [10,  11].  Wherever  prominent  three-fold  symmetrical  colony  structures  of 
c-Zr02  were  observable  on  transverse  sections,  there  were  notable  angular  differences  between 
the  symmetry  axes  of  these  colonies  indicating  the  presence  of  small  angle  tilt  boundaries  in 
the  AI2O3,  (often  in  the  range  of  10°)  with  tilt  axes  parallel  to  the  growth  axis.  Since  such  tilt 
boundaries  would  not  be  stressed  in  samples  under  tension  parallel  to  the  bar  axis,  they  were 
not  considered  to  be  of  any  importance  in  the  creep  behavior. 

The  c-Zr02  phase,  in  turn,  exhibited  various  morphological  shapes  ranging  from  aligned 
fibers  or  platelets  in  either  well  organized  colonies  or  in  less  well  ordered  forms,  usually  had 

1  About  4.2  mol%  of  Y2O3  is  required  to  form  cubic  zirconia.  This  was  the  case  with  all  the  eutectics  in  the 
present  study.  Therefore,  hereafter  the  zirconia  is  referred  to  only  as  c-ZrC>2 
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came  from  the  Mechanical  Engineering  Department  at  M.I.T.  in  some  equipment  acquisitions 
and  from  a  NATO  Fellowship  to  a  visiting  scientist. 
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one  narrow  dimension  in  the  0.2/xm  range,  and  was  always  fully  encapsulated  in  the  AI2O3 
phase  with  high  quality  nearly  coherent  interfaces  [13],  also  verified  in  the  present  work  with 
several  direct-lattice-imaging  TEM  observations  [23].  The  eutectic  had  fairly  reproducible  and 
reasonably  attractive  levels  of  fracture  toughness  [12]. 

Corresponding  determination  of  the  texture  of  the  c-ZrC>2  phase,  largely  based  on  electron 
diffraction  information,  indicated  that  this  phase  had  primarily  a  <  112  >  growth  texture 
with  fiber  symmetry,  incorporating  random  rotations  about  this  axis1  Figure  1  shows  a  typical 
morphology  of  the  A^Os/c-ZrC^  eutectic.  The  sketch  of  Fig.  2  shows  pictorially  the  mutual 
orientational  association  of  the  two  component  phases. 

The  DS  eutectics  contained  significant  levels  of  residual  stress  at  room  temperature,  resulting 
from  thermal  misfit  between  the  component  phases.  The  AI2O3  was  in  a  state  of  uniaxial  tension 
parallel  to  the  axis,  at  a  level  of  1.0  GPa  while  the  c-ZrC>2  was  in  some  form  of  multiaxial 
compression,  roughly  at  the  same  level  [10].  Such  residual  stresses  that  decrease  in  level  at 
elevated  temperatures  were  rapidly  relaxed  on  the  onset  of  steady  creep  flow. 

2.2  Creep  equipment  and  creep  strain  measurement 

The  creep  experiments  were  carried  out  in  tension  inside  a  Centorr  vacuum  chamber  containing 
a  specially  designed  hot  zone  of  5cm  diameter  and  5cm  length  made  up  of  a  split  cylindrical 
configuration  of  a  1mm  thick  Ta  sheet  with  alternating  up-and-down  slits  to  increase  the  path 
length  of  the  heating  element.  The  hot  zone  was  surrounded  by  a  series  of  Mo  radiation  shields. 
The  long  specimens  threaded  through  the  hot  zone  and  were  gripped  at  their  ends  by  specially 
produced  stainless  steel  friction  grips,  clamping  down  on  the  ends  of  the  specimens  through 
annealed  soft  nickel  sheets  to  avoid  high  local  contact  forces  on  the  rods  that  could  have  resulted 
in  fracture.  Two  relatively  massive  water  cooled  OFHC  Cu  plates  were  placed  between  the 
two  exits  of  the  hot  zone  and  the  specimen  grips  to  assure  that  the  gripping  conditions  were 
not  to  exceed  about  200°C  to  avoid  slippage  in  the  grips.  Figure  3  shows  a  schematic  layout  of 
the  interior  of  the  high  temperature  high  vacuum  creep  chamber. 

An  optical  grade  sapphire  window  permitted  viewing  the  hot  specimen  through  a  narrow 
axial  slit  in  the  heating  elements  to  permit  direct  measurements  of  the  specimen  temperature 
by  means  of  a  two-wave-length  pyrometer  (from  Omega  Vanzetti,  Sharon,  MA.),  providing 
emissivity-independent  measurements.  In  addition  to  a  control  thermocouple  inside  the  hot 
zone,  another  thermocouple  placed  close  to  the  specimen  without  touching  it,  was  used  to 
actually  record  the  steady  temperature  in  the  thermal  cavity.  The  temperature  measured  by 
this  second  thermocouple  and  that  measured  by  the  pyrometer  on  the  specimen  usually  agreed 
quite  well.  The  hot  zone  temperature  was  automatically  controlled  by  the  PID  control  system 
of  the  Centorr  equipment. 

1  Textures  such  as  <  111  >  have  also  been  reported  by  other  investigators. 
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Since  the  steady  state  creep  extensions  of  the  specimens  under  stress  were  generally  quite 
small,  to  eliminate  random  flexures  and  motions  of  the  massive  creep  frame  from  influencing 
the  measurements,  the  displacements  of  the  two  ends  of  the  load  train  immediately  outside  the 
vacuum  chamber  were  simultaneously  measured  by  sets  of  four  LVDTs  at  both  the  top  and 
bottom  ends  of  the  Centorr  chamber.  The  differences  between  these  two  measurements  then 
corresponded  to  the  extension  of  the  specimen.  Since  the  temperature  of  the  specimen  inside 
the  hot  zone  was  higher  by  1200-1300°C  above  the  portions  of  the  sample  outside  the  hot  zone, 
the  gauge  length  was  taken  as  the  axial  extent  of  the  hot  zone.  Even  under  ideal  conditions  the 
total  creep  strain  was  never  too  uniquely  determinable.  To  overcome  this  difficulty  all  creep 
strains  at  steady  state  were  always  measured  incrementally  at  a  given  temperature  under  two 
different  applied  stresses  and  often  as  loading  and  unloading  cycles.  This  practice  of  relative 
measurements  of  creep  strain  did  not  only  give  more  reliable  determinations  of  the  strain  rates  at 
the  two  different  levels  of  stress  but  also  demonstrated  the  nearly  complete  absence  of  transients 
in  such  incremental  changes  in  steady  state  creep.  A  typical  response  of  this  type  is  shown  in 
Fig.  4  for  creep  at  1520°C  at  stress  levels  of  150MPa  and  200MPa.  The  accuracy  of  strain 
measurement  was  around  1.4  x  10~5. 

2.3  Fracture  toughness  experiments 

Fracture  toughness  of  the  eutectics  was  measured  by  the  well  known  technique  of  determining 
the  extent  of  cracking  from  the  corners  of  Vickers  hardness  indentations,  but  by  utilizing  a 
more  exact  technique  of  determining  the  open  crack  contours  at  the  surfaces  as  described  by 
Haubensack  and  Argon  [14].  The  fracture  toughness  was  also  estimated  from  the  size  and  shape 
of  pores  in  the  as-grown  eutectics  when  they  gave  rise  to  premature  fracture.  Some  estimates 
were  also  made  from  the  measured  toughness  profiles  of  fracture  surfaces. 

2.4  Computational  methodologies 

Two  separate  computational  studies  were  carried  out  in  the  program.  In  the  one  carried  out  on 
the  dislocation  core  structures  in  sapphire  at  the  atomic  level,  the  MOLDY  MD  programs  were 
used  to  calculate  core  structures  of  edge  dislocations  on  basal  planes  and  pyramidal  planes  of 
AI2O3.  In  the  FEM  computations  of  distributions  of  stresses  among  the  two  different  compo¬ 
nents  of  AI2O3  and  c-Zr02  during  creep,  the  commercial  ABAQUS  code  was  used. 
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3  Experimental  Results 

3.1  Creep  experiments 

3.1.1  Transient  effects 

Upon  initial  application  of  stress  to  a  virgin  specimen  a  strain  transient  was  always  observed. 
One  such  transient  recorded  for  an  experiment  at  1200°C  under  a  stress  of  300  MPa  is  shown  in 
Fig.  5.  The  characteristic  time  constant  of  the  transient  decreased  with  increasing  temperature, 
e.g.  from  c.a.  3.4  x  103sec  at  1200°C  to  200  sec  at  1400°C.  The  usual  explanation  of  the 
transient  to  be  due  to  initial  work  hardening  in  the  creeping  components  prior  to  establishment 
of  a  steady  state  by  recovery  processes,  familiar  in  creep  of  homogeneous  metals,  was  discounted 
here  since  no  prominent  transients  were  found  following  stress  increases  at  steady  state  creep. 
The  most  likely  cause  of  the  initial  transient  was  the  rapid  stress  relaxation  in  the  coarse  c-ZrC>2 
component  by  creep  since  for  the  measured  texture  of  that  component  there  is  an  abundance  of 
the  <  110  >  {100}  principal  slip  systems  that  are  well  oriented  for  slip,  and  the  creep  resistance 
of  c-ZrC>2  in  this  temperature  range  is  known  to  be  poor  [15].  In  more  ideal  morphologies,  as 
presented  by  Sayir  and  Farmer  [2]  and  Argon,  et  al.  [16],  the  coarse  fraction  of  the  c-Zr02  is 
located  on  colony  borders.  In  the  usual  less  ideal  morphologies  where  fibrillar  c-ZrC>2  is  still 
aggregated  but  does  not  form  well  structured  colonies  the  coarse  fraction  was  found  to  be  less 
regularly  distributed  (see  Fig.  1).  An  upper  bound  manifestation  of  the  poor  creep  resistance 
of  c-ZrC>2  can  be  assessed  from  a  Finite  Element  (FEM)  study  that  will  be  presented  in  Section 
4.1.  There  it  was  found  that  for  the  complete  relaxation  of  the  deviatoric  stresses  in  the  entire 
c-ZrC>2  phase  the  volume  average  deviatoric  stresses  as  well  as  mean  normal  stresses  in  the 
AI2O3,  increase  roughly  by  25-30%  which  would  result  in  an  increment  of  additional  elastic 
strain  in  the  sample  of  about  6  x  10-4.  This  level  is  illustrated  in  Fig.  5  by  the  horizontal  line. 
The  actual  amount  of  stress  relaxation  in  the  c-ZrC>2,  however,  is  difficult  to  determine  since  a 
substantial  fraction  of  this  component  has  phase  dimensions  in  the  range  of  0.2-0.4//m  and  is 
likely  to  be  dislocation-free  and  incapable  of  plastically  deforming  as  we  discuss  in  Section  5.0. 
These,  and  the  difficulties  in  determining  the  absolute  measures  of  strain  derived  from  the  hot 
zone  portions  of  the  samples  requires  us  to  de-emphasize  this  portion  of  the  creep  response. 

3.1.2  Steady  state  creep 

As  already  mentioned  in  Section  2.2  the  main  information  on  steady  state  creep  was  obtained 
from  incremental  experiments  of  sudden  stress  increases  and  decreases  at  constant  temperature 
as  shown  in  the  typical  case  of  Fig.  4  of  creep  response.  A  simple  check  of  the  magnitude 
of  the  instantaneous  stretches  or  contractions  for  stress  increases  or  decreases,  utilizing  the 
appropriate  information  on  the  temperature  dependent  Young’s  modulus  of  AI2O3  [17]  gave 
in  all  cases  that  the  recorded  jumps  were  all  about  a  factor  of  2  larger  than  what  could  be 
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expected  from  the  sample  in  the  hot  zone.  Since  this  was  well  within  the  additional  flexures 
and  relaxations  in  the  load  train  the  recorded  strain  rates  were  considered  to  be  reliable. 

After  the  final  configuration  of  the  hot  zone  discussed  in  Section  2.2  was  completed,  11 
determinations  of  steady  state  creep  rate  were  made,  all  together  at  3  temperatures  of  1200, 
1400,  1520°C  and  at  stress  levels  ranging  from  150MPa  to  275MPa,  derived  from  incremental 
experiments  similar  to  those  in  Fig.  4.  These  measurements  are  listed  in  Table  1  and  are 
plotted  in  Fig.  6  together  with  the  creep  rates  at  1400°C  due  to  Sayir  and  Farmer  [2],  also 
quoted  in  [16] 1.  The  stress  exponents  of  the  creep  rates  at  1400°C  and  1520°C,  shown  in  Fig. 
6,  are  4.16  and  5.38  and  are  lower  than  that  of  Sayir  at  6.00.  As  we  will  discuss  in  the  creep 
model  in  Section  5.1,  these  exponents  are  relatively  high  and  need  to  be  explained,  as  we  will 
do  later  in  Section  5.2. 

Table  1  contains  also  sufficient  information  for  the  determination  of  the  activation  energy 
of  the  governing  creep  process.  The  plot  of  the  creep  rates  at  a  stress  of  200MPa  for  three 
temperatures  of  1200°C,  1400°C  and  1520°C  in  Fig.  7  gives  an  activation  energy  of  Q  = 
71.1  kcal/mol  for  the  rate  controlling  process  in  steady  state  creep  which  we  consider  to  be 
due  to  oxygen  ion  diffusion  that  is  known  to  be  the  slowest  diffusing  species  in  AI2O3.  This 
value  is  to  be  compared  with  105  kcal/mol  for  temperatures  above  1400°C  and  26  kcal/mol  for 
temperatures  below  1400°C  for  the  same  process  measured  by  Oishi  and  Kingery  [18]  in  pure 
polycrystalline  AI2O3.  The  somewhat  lower  magnitude  of  the  measured  activation  energy  for 
oxygen  ion  diffusion  is  attributed  to  the  likelihood  that  a  significant  fraction  of  the  diffusion 
current  must  have  been  along  the  the  incoherent  interfaces. 


3.2  Fracture  toughness  measurements 


As  mentioned  in  Section  2.3,  the  fracture  toughness  was  estimated  by  the  open  crack  contours 
at  the  Vickers  indentation  corners.  Fig  8  shows  the  Vickers  indentation  and  the  open  crack 
extending  from  one  of  the  corners.  According  to  Haubensack  and  Argon  [14],  the  K jc  can  be 
estimated  from 

{COD)E 


KIC  = 


F(r/a)r1/2 


(1) 


where  F(r/a)  is  a  non-dimensional  characteristic  crack-shape  factor  dependent  on  the  geometry 
and  the  mode  of  loading,  a  is  the  crack  length,  and  E  is  the  Young’s  modulus.  The  estimated 
values  of  Kic  of  eutectic  Al203/c-Zr02  ranged  from  2.0  MPam1/2  to  5.0  MPam1/2. 

The  fracture  toughness  can  also  be  estimated  from  the  size  and  shape  of  pores  in  the  as- 
grown  eutectic  and  surface  flaws  when  they  lead  to  premature  fracture.  In  Fig  9,  an  internal 
pore  and  a  surface  flaw  are  identified.  The  estimated  fracture  toughness  is  around  1.0  MPam1/2. 


1The  higher  creep  rates  and  the  larger  stress  exponents  in  these  early  investigations  are  most  likely  a  result 
of  not  reaching  steady  state  but  still  being  in  a  late  stage  of  transient  behavior. 
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The  crack  deflection  by  the  morphology  can  be  a  very  effective  way  in  toughening.  By 
measuring  the  fracture  surface  profile,  the  energy  release  rate  can  be  estimated.  In  Fig  10  the 
AFM  surface  profile  characterization  of  a  fractograph  is  shown.  It  suggests  that  the  critical 
energy  release  rate  Gico  is  around  1.0  J/m2,  which  is  reasonable. 


4  Computational  Results 

4.1  Distribution  of  stress  in  the  Al203/c-ZrC>2  eutectics  in  creep 

A  very  important  way  in  which  creep  in  the  eutectics  differs  from  creep  in  single  crystal  sapphire 
is  that  in  the  former  the  stresses  are  distributed  in  a  complex  manner.  In  addition  to  resid¬ 
ual  stresses  arising  from  the  different  thermal  expansions  of  the  two  components,  the  applied 
stresses  result  in  complex  internal  local  stress  distributions  due  to  different  elastic  properties  of 
the  A1203  and  the  c-Zr02  component  phases.  Thus,  to  develop  some  necessary  understanding 
of  this  phenomenon  on  a  broad  basis,  a  linear-elastic  FEM  analysis  was  carried  out  in  a  2-D 
plane  strain  setting  of  the  eutectic  in  the  ordered  regions  shown  in  Fig.  1.  Here  the  c-Zr02 
domains  in  2-D  are  idealized  as  elliptical  cylinders  occupying  a  volume  fraction  of  0.33  as  de¬ 
picted  in  Fig  11.  The  chosen  representative  computational  volume  (RVE)  is  indicated  in  this 
figure.  We  have  considered  the  c-Zr02  to  be  always  fully  relaxed  of  all  shear  stresses  but  with 
an  unrelaxed  appropriate  bulk  modulus  of  83GPa.  This  is  accomplished  operationally  in  the 
FEM  code  for  purely  elastic  behavior  by  choosing  the  Young’s  modulus  of  the  c-Zr02  artifically 
as  0.05MPa  and  the  Poisson’s  ratio  as  0.4999999.  All  other  chosen  material  constants  are  listed 
in  Table  2.  The  broader  analysis  was  considered  for  several  aspect  ratios  a/b  of  ellipsoidal  rods 
of  c-Zr02  as  1,3,4  and  5.  Of  these  results  we  present  here  only  the  distributions  for  a/b  =  3.0 
as  most  representative.  All  analyses  were  limited  to  the  elastic  range.  In  this  FEM  framework 
we  considered  several  process  simulations:  a)  an  applied  tensile  stress  az  alone;  b)  cooling 
the  morphology  from  a  stress-free  condition  at  1875°C  down  to  1400°C  to  determine  residual 
stresses  at  1400°  and;  c)  cooling  as  in  (b)  plus  an  applied  tensile  stress  of  az  =  300MPa.  Of 
these,  we  discuss  here  only  the  case  of  the  internal  stress  distribution  under  an  applied  stress 
crz  =  300MPa,  and  the  result  of  a  cooling  history,  from  1875°C  (Since  behavior  is  linear-elastic, 
results  for  all  other  applied  stress  levels  can  be  determined  by  re-scaling). 

In  Fig.  12a  we  show  the  distribution  of  the  Mises  de viator ic  stresses  for  the  case  of  az  = 
300MPa  and  T=1400°C,  in  both  the  A1203  and  c-Zr02  regions.  As  expected,  the  states  of 
stress  inside  the  ellipsoidal  regions  of  the  c-Zr02  are  close  to  constant.  Figure  12b  gives  the 
more  important  distribution  of  the  climb  stress  aj ,  parallel  to  the  pyramidal  dislocation  Burgers 
vector  at  an  angle  of  <j>  =  57.7°  as  shown  in  Fig.  17.  The  volume-average  level  of  this  important 
climb  producing  stress  was  found  to  be  0\  =  216.3MPa  which  gives  the  factor  q=  o\/az  to  be 
0.721.  This  value  appears  surprisingly  low  in  view  of  the  assumed  complete  relaxation  of  the 
deviatoric  stresses  in  the  c-Zr02,  until  it  is  recognized  that  the  c-Zr02  still  supports  fully  a 
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mean  normal  stress  which  limits  the  stress  enhancement  in  the  AI2O3  component.  The  factor 
q  enters  in  an  important  way  into  the  creep  model  of  the  eutectics  as  described  in  Section  5.2. 

Figure  13a  and  13b  show  the  corresponding  results  for  internal  deviatoric  (von  Mises) 
stresses  and  mean  normal  stresses  at  1400°C,  resulting  from  a  cooling  history  from  1875°C 
to  1400°C.  Examination  of  Figs.  13a  and  13b  shows  that  the  thermal  misfit  in  the  AI2O3  at 
1400°C  are  in  the  range  of  substantial  fractions  of  one  GPa,  and  clearly  consistent  with  the 
reported  residual  stresses  at  room  temperature  referred  to  in  Section  2.1. 

4.2  Simulations  of  the  core  structures  of  edge  dislocations  in  AI2O3 

As  is  expanded  further  in  the  creep  model  presented  in  Section  5.1  the  creep  resistance  of  the 
eutectics  are  governed  by  that  of  the  topologically  continuous  AI2O3  component  which  has  a 
remarkably  tight  growth  texture  of  [0001]  parallel  to  the  growth  axis  of  the  eutectic  rods.  This 
renders  the  best  glide  systems  in  the  base  plane  and  even  the  prismatic  system  inoperative. 
Thus,  for  geometrical  reasons  extension  of  the  AI2O3  in  the  [0001]  direction  by  glide  or  climb 
can  only  be  possible  through  the  motion  of  the  pyramidal  dislocations.  According  to  general 
consensus  glide  deformation  on  the  pyramidal  dislocations  is  too  difficult  to  happen  and  has 
never  been  reported,  leaving  climb  of  such  dislocations  as  the  only  realistic  possibility.  To 
establish  the  basis  of  the  impossibility  of  glide  on  the  pyramidal  system  a  molecular  dynamics 
(MD)  simulation  was  undertaken  of  the  core  structure  of  edge  dislocations  1/3  <  1101  > 
on  this  slip  system  using  the  established  MOLDY  [26]  code.  For  comparison  purposes  the 
core  structure  of  the  basal  dislocations  1/3  <  2110  >  was  also  simulated.  The  full  details  of 
the  simulation  results  which  are  beyond  the  scope  of  this  final  report  have  been  presented  in 
a  comprehensive  communication  [5]  and  will  not  given  here.  Here  only  the  most  important 
findings  are  presented. 

Figure  14  shows  the  cation  sub-lattice  of  AI2O3  which  establishes  the  principal  dislocation 
Burgers  vectors  based  on  the  periodicity  of  the  A1  cation  holes.  The  three  principal  slip  systems: 
in  the  base  plane;  in  the  prismatic  system  and  finally  in  the  pyramidal  system  are  clearly  iden¬ 
tified.  The  simulations  of  the  base  plane  dislocations  and  the  pyramidal  dislocations  required 
selection  of  two  separate  simulation  cells  shown  in  Fig.  15  which  gives  the  crystallographic 
symmetry  of  the  idealized  AI2O3  lattice  as  viewed  from  the  ”c”  direction.  In  the  figure  the 
large  open  circles  show  the  hexagonal  packing  of  the  Oxygen  anions  in  the  base  planes  whereas 
the  small  solid  black  circles  and  the  small  open  circles  show  the  location  of  the  A1  cations  and 
cation  holes  respectively,  on  an  interstitial  plane  immediately  above  the  Oxygen  plane. 

The  rectangle  ABCD  gives  the  geometry  of  the  base  of  the  simulation  cell  chosen  to  simulate 
the  cores  of  the  1/3  <  2110  >  basal  dislocation,  while  the  rectangle  A’ACC’  gives  the  geometry 
of  the  base  of  the  simulation  cell  chosen  for  simulation  of  the  cores  of  the  1/3  <  1101  > 
pyramidal  dislocation.  In  both  cases  periodic  boundary  conditions  were  used  to  obtain  valid 
simulations  of  edge  dislocation  cores,  unaffected  by  surface  conditions.  To  obtain  realistic 
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conditions  with  rapid  convergence  of  the  MD  energy  relaxations  real  atom  coordinates  for  the 
ion  lattices  were  used  which  differ  substantially  from  the  ideal  coordinates  presented  initially 
by  Kronberg  [27].  In  both  simulations  to  avoid  unwanted  long  range  lattice  curvatures  edge 
dislocation  dipoles  were  introduced  into  the  simulation  cells  according  to  a  rigid  prescription 
discussed  elsewhere  [5]. 

The  principal  finding  related  to  the  immobility  of  the  pyramidal  dislocation  is  shown  in 
Fig.  16a  and  16b  which  views  the  AI2O3  lattice  from  the  [1120]  direction  (refer  to  Fig.  14). 
Figure  16a  shows  the  unrelaxed  core  of  the  initial  configuration  of  the  1/3  <  1101  >  pyramidal 
dislocation  located  in  the  (1102)  planes  made  up  entirely  of  cation  holes  in  the  cation  sub¬ 
lattice.  The  sketch  on  the  right  side  of  the  figure  delineates  the  two  extra  ion  half  planes  of 
the  dislocation.  The  final  relaxed  configuration  of  the  core  after  minimization  of  system  energy 
by  the  MOLDY  program  is  shown  in  Fig.  16b.  The  figure  indicates  that  the  form  of  energy 
minimization  that  occurs  and  affects  the  core  structure  splits  the  core  into  two  fractional  edge 
dislocations  lying  on  parallel  cation  planes  on  either  side  of  the  initially  more  favorable  plane 
of  cation  holes,  being  displaced  to  planes  with  both  high  level  covalent  and  ionic  bonding  as 
delineated  in  the  sketch  on  the  right  side  of  the  figure.  Thus,  while  the  core  energy  is  distinctly 
reduced  the  two  partial  edge  dislocations  are  rendered  sessile  as  a  detailed,  so-called,  gamma 
surface  topography  of  the  potential  shear  resistance  across  these  planes  demonstrate. 

The  corresponding  simulation  of  the  core  structure  of  the  base  plane  edge  dislocation  1/3  < 
2110  >  which  is  not  presented  here  exhibits  a  very  similar  splitting  of  the  core  structure  of  the 
dislocation  as  in  the  case  of  the  pyramidal  dislocation.  In  this  case,  however,  since  each  base 
plane  between  adjacent  close  packed  anion  planes  is  identical  to  all  other  planes  no  reduction 
of  mobility  results  from  the  core  rearrangement. 

Thus,  in  summary,  the  MD  simulations  have  clearly  demonstrated  that  the  pyramidal  dis¬ 
location  can  not  glide  but  that  their  climb  motion  should  be  relatively  unimpaired. 

5  The  Creep  Model 

5.1  Basic  assumptions 

The  fact  that  the  c-ZrC>2  of  considerably  lower  creep  resistance  is  encapsulated  everywhere  by 
the  topologically  continuous  AI2O3  with  a  very  tight  growth  texture  of  [0001]  along  the  axes  of 
the  DS  eutectic  rods  indicates  that  steady  state  creep  must  be  controlled  by  the  AI2O3  phase, 
once  any  limited  compliance  increment  due  to  the  stress  relaxation  in  c-ZrC>2  is  complete.  This 
suggests  that  the  overall  creep  response  must  have  the  same  characteristics  of  steady  state 
creep  in  sapphire  single  crystals  as  reported  by  Firestone  and  Heuer  [3],  albeit  with  certain 
important  added  complexities  derived  from  the  presence  of  the  c-ZrC>2  phase.  Figure  1  shows 
the  well  known  phase  morphology  of  a  typical  Al203/c-Zr02  eutectic  where  the  bright  phase 
is  the  c-ZrC>2.  Very  much  like  the  7’  phase  in  superalloy  single  crystals  of  CMSX-3  [18],  a 
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large  fraction  of  the  c-Zr02  has  dimensions  in  the  range  of  0.2-0.4  //m  and  must  be  regarded 
as  being  too  small  to  have  the  capability  of  undergoing  plastic  flow  by  independent  internal 
processes  of  dislocation  multiplication  by  glide  or  climb.  The  c-Zr02  phase  components  of 
larger  dimensions  (the  irregular  shaped  ones  in  Fig.  1)  should  undergo  rapid  relaxation  of 
shear-  stresses  in  the  1200-1520°C  temperature  range  of  interest,  as  the  transient  creep  curve 
of  Fig.  5  suggests.  However,  as  the  FEM  analysis  of  Section  4.1  shows,  complete  shear  stress 
relaxation  in  the  c-Zr02  still  leaves  the  mean  normal  stress  intact  in  it  which  prevents  large 
average  stress  enhancement  in  the  A12C>3  phase. 

Since  the  best  slip  systems  of  the  basal  and  prismatic  type  are  largely  unstressed  for  glide  due 
to  the  tight  texture  of  the  Al203  phase,  and  since  the  pyramidal  dislocations  of  1/3  <  1101  > 
type  are  sessile  in  glide,  as  discussed  in  Section  4.2,  creep  in  this  component  for  eutectic  bars  in 
tension,  can  only  be  a  consequence  of  climb  of  the  (1/3)  <  1101  >  pyramidal  edge  dislocations. 
Nabarro  [19]  has  presented  an  idealized  model  of  creep  derived  entirely  from  climb  of  dislocations 
of  a  Frank-type  network  where  material  fluxes  are  between  the  actual  climbing  dislocations  with 
Burgers  vectors  parallel  to  the  stress  axis  and  those  of  other  types  experiencing  no  climb  forces. 
Firestone  and  Heuer  [3]  found  reasonable  agreement  of  their  results  with  the  Nabarro  model 
where  the  stress  exponents  were  in  the  range  of  3  as  that  model  predicts.  In  the  case  of  the 
eutectics  where  the  stress  exponents  are  in  the  range  of  4.5-6. 0  different  conditions  should  hold, 
requiring  considerable  modification  of  the  Nabarro  framework. 

In  our  creep  model  it  is  assumed  that  the  (1/3)  <  1101  >  pyramidal  dislocations  pre-exist 
in  the  A12C>3  phase,  and  will  multiply  by  well  known  topological  convolution  processes  from 
Bardeen-Herring  [20]  sources,  and  will  climb  under  stress  during  creep.  The  mean-free-path 
lengths  of  climbing  dislocations  in  Al203  should  involve  several  multiples  of  the  interphase 
dimensions,  requiring  them  to  repeatedly  bow  and  straighten-out  as  they  thread  through  the 
small  c-Zr02  domains  that  will  remain  impenetrable  to  them.  Such  repetition  of  transient 
line  shapes  can  be  viewed  as  the  primary  source  of  the  elevation  of  the  creep  stress  exponents 
to  become  larger  than  3  [21].  The  climb  of  these  dislocations,  will  be  diffusion  controlled, 
provided  that  adequate  jog  concentrations  are  present  along  the  dislocations.  Such  jogs  are 
expected  to  be  produced  by  nodal  emission  processes  from  network  junctions  rather  than  being 
of  a  thermal  equilibrium  nature  which  preliminary  analyses  show  to  be  far  too  difficult.  The 
diffusional  material  transport  will  be  between  dislocation  cores  and  the  large  concentrations 
of  incoherent  or  semi-coherent  interfaces.  These  interfaces  will  be  completely  opaque  to  the 
transmission  of  dislocations  from  c-Zr02  to  A1203  or  vice  versa.  Residual  stresses  due  to  initial 
thermal  misfit  between  the  components  will  be  relieved  during  the  transient  phases  of  the  creep 
and  provide  another  contribution  to  the  recorded  creep  strain  shown  in  Fig.  5,  but  that  steady 
state  creep  will  be  governed  by  the  climb  forces  resulting  entirely  from  the  applied  stress. 
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5.2  The  creep  rate 


Consider  a  round  tensile  creep  bar  of  an  Al203/c-Zr02  eutectic  as  sketched  in  Fig.  17  where 
local  axes  1,2,3,  are  chosen  in  a  pyramidal  glide  system  with  axis  1  parallel  to  the  1/3  <  1101  > 
Burgers  vector,  axis  2  normal  to  the  pyramidal  plane  making  an  angle  with  the  bar  axis,  and 
axis  3  parallel  to  the  line  vector  of  a  positive  edge  dislocation.  Under  a  uniform  stress  <j\  a 
positive  edge  dislocation  will  climb  in  the  negative  2  direction  with  a  velocity  vc  to  contribute  to 
a  uniaxial  tensile  creep  strain  rate  parallel  to  the  1  direction.  This  velocity,  if  it  were  governed 
by  diffusion  controlled  climb  would  be  given  by  [22]: 


27rD(<7i  —  (TT)b 2 
kTln(rs/rc ) 


(2) 


where  D  =  D0exp(—Q / RT)  is  the  diffusion  constant  of  AI2O3  where  oxygen  is  recognized  to 
be  the  slow  diffusing  species,  ot  is  the  triaxial  component  of  the  applied  stress  (=  crz/3)  that 
does  not  promote  equilibrium  climb,  rs  is  the  distance  from  the  dislocation  core  to  vacancy 
sinks  along  the  interfaces,  rc  the  core  radius  from  which  point  defect  emission  is  considered,  b 
the  magnitude  of  the  Burgers  vector,  £>0  the  pre-exponential  factor  of  the  diffusion  constant, 
Q  the  activation  energy  of  O  ion  diffusion  in  AI2O3  and  R  and  T  have  their  usual  meaning. 
The  kinematics  of  shape  change  due  to  climb  alone  of  dislocations  has  been  considered  more 
broadly  by  Groves  and  Kelly  [23].  In  the  case  of  the  eutectics  a  more  specialized  consideration 
is  adequate.  Thus,  considering  that  there  will  be  three  equivalent  pyramidal  systems  that  can 
contribute  to  the  axial  strain  rate  equally,  the  axial  creep  strain  rate  ez  along  the  bar  and  the 
accompanying  radial  contractile  strain  rates  er  are  given  respectively  as  [24]: 


and 


where 


ez  =  (2  —  3  cos20)eo 
3 

er  =  —  (1  —  -cos2<j))e 0 


6q  —  bpmvc 


(3a) 

(3b) 

(4) 


is  the  main  creep  rate  in  the  pyramidal  system  coordinates  without  regard  to  overall  volume 
preservation  ,  which  is  considered  to  result  in  Eqns(3a)  and  (3b)  [24].  Relating  the  actual  climb- 
producing  tensile  stress  <7i  to  the  axial  stress  <rz  and  considering  the  volume-average  effect  of 
the  local  variation  of  stress  due  to  the  presence  of  the  c-Zr02  components  by  a  factor  q  gives 


<Ji-crT  =  ^qaz(  1  -  cos2<f>) 


(5) 
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If  the  mobile  dislocation  density,  pm,  of  climbing  dislocations  at  steady  state  is  governed  by 
mutual  interactions  in  a  self  adjusting  basis  [19];  i.e. 

Pm  =  (^)V(l  -  cosW  (6) 

one  has,  by  finally  combining  Eqns  (3a), (4), (5)  and  (6)  the  axial  steady  state  creep  rate: 

_  \/0z\3(PQ\  (1  ~  cosV)3(2  —  3 cos  </>) 

z~A[b2’[  n}  [kT}  ln(r3/rc)  U 

where 

A  =  (8a) 

and 

/3  =  n/b 3  =  0.079  (8b) 

where  is  the  ionic  volume  of  O  in  the  AI2O3  lattice  and  q  is  a  factor  which  relates  the  local 

volume-average  climb  stress  to  the  axial  stress  az,  as  was  determined  from  the  FEM  analysis 
of  the  stress  distribution  in  the  AI2O3  phase  under  an  applied  tensile  stress,  as  presented  in 
Section  4.1. 

Eqn.  (7)  gives  the  creep  rate  due  to  the  climb  of  quasi-straight  dislocations  in  a  homogeneous 
stress  field.  As  is  discussed  in  Section  5.3  this  is  not  the  case  in  the  Al203/c-Zr02  eutectics 
where  the  climbing  dislocations  need  to  thread  through  the  non-deforming  c-ZrC>2  domains 
acting  as  dispersoids  and  are  required  to  alternately  bow  around  these  domains  and  be  released 
to  straighten  out,  acting  effectively  as  if  straight  dislocations  were  moving  through  a  strongly 
varying  internal  stress  field  which  will  increase  the  stress  exponent  in  predictable  ways  and 
decrease  the  actual  creep  rate  [21]. 

In  Section  5.4  the  creep  model  is  compared  with  the  experimental  results. 

5.3  Effect  of  internal  resistance  variations 

As  stated  in  Section  5.2  the  creep  response  of  the  Al20s/c-Zr02  eutectics  differ  significantly 
from  the  creep  behavior  in  homogeneous  sapphire  single  crystals  studied  by  Firestone  and  Heuer 
[3].  First  because  of  the  non-uniform  stress  distribution  in  the  morphology  of  AI2O3  and  c-ZrC>2 
discussed  in  Section  4. 1  above  but  even  more  importantly  by  the  large  perturbations  that  the 
climbing  dislocations  encounter  threading  through  the  isolated  and  largely  non-deforming  c- 
ZrC>2  domains.  Thus,  consider  the  convolutions  that  a  climbing  dislocation  has  to  go  through 
as  it  threads  through  the  gaps  of  the  c-ZrC>2  domains  as  depicted  in  Fig.  18.  First,  the  climbing 
dislocation  must  squeeze  into  the  gaps  between  the  domains  as  depicted  in  Fig.  18a  to  a  critical 
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configuration  much  like  the  Orowan  bowing  process  for  non-shearable  dispersoids.  The  peak 
stress  that  is  required  for  this  configuration  to  be  achieved  is 


a  =  a i  — 


2£ 

bA 


fj,b 

1 


(9) 


where  £  is  the  dislocation  line  tension  (=  /i62/2)  and  A  is  the  size  of  the  inter-domain  gap. 
Clearly,  here  the  process  is  not  one  of  glide  but  climb  where  a  is  the  climb  stress.  Once 
the  critical  configuration  is  reached  and  the  climbing  dislocation  surrounds  the  impenetrable 
domain  and  pinches  off  as  shown  in  Fig.  18b,  the  separated  dislocation  has  acquired  a  definite 
cusp  where  this  shape  of  the  dislocation  will  undergo  accelerated  climb  due  to  the  advantageous 
line  tension  effects.  These  required  contortions  of  the  dislocations  are  viewed  as  if  they  were 
the  same  as  positive  and  negative  internal  stresses  cr*,  retarding  and  then  accelerating  climb, 
very  much  like  the  corresponding  cases  of  favorable  and  unfavorable  internal  shear  stresses  that 
a  glide  dislocation  needs  to  go  through  between  dispersoids.  Such  problems  were  treated  in 
great  detail  by  many  investigators,  but  particularly  thoroughly  by  Li  [21],  which  illustrated 
that  a  gliding  dislocation  (in  the  present  case  a  climbing  dislocation)  loses  more  time  in  the 
region  of  adverse  internal  stress  than  the  time  it  gains  in  moving  through  regions  of  favorable 
internal  stress.  The  effect  is  illustrated  in  Fig.  19  where  19a  shows  the  sinusoidally  varying 
internal  stress  <jj  and  the  prevailing  climb  stress  a.  The  most  important  consequence  of  such 
deceleration  and  acceleration  of  motion  is  an  increase  in  the  effective  stress  exponent  m  of  the 
dislocation  velocity,  scaled  by  the  ratio  of  a  i/a  shown  in  Fig.  19b.  A  secondary  effect  is  a  factor 
C  =  v/v  that  gives  the  decrease  in  the  average  dislocation  velocity  relative  to  the  velocity  in 
unhindered  climb  shown  in  Fig.  19c,  also  scaled  with  the  ratio  of  (Tj/cr.  To  assess  a  measure 
of  the  effect  it  must  be  noted  from  Fig.  1  that  in  much  of  the  morphology  the  inter-domain 
distance  A  would  appear  to  be  between  0.3-1.0/im.  Taking  the  small  dimension  to  determine 
the  maximum  level  of  this  effect,  together  with  b  =  5.12  x  10-10m  and  fx  —  E/2{\+v)  =  73 GPa 
[17]  at  1400°C  we  determine  a*  =  125MPa.  For  an  applied  stress  oz  =  300MPa,  in  the  range 
of  interest,  that  would  give  a  volume-average  climb  stress  a  =  qaz(  1  —  cos2<f>)  of  153MPa,  and 
a  ratio  «  0.815,  Fig  19b  indicates  that  the  effective  exponent  of  the  stress  in  the  velocity 
expression  should  increase  to  nearly  m  =  3  from  unity  as  the  maximum  effect.  Moreover,  Fig. 
19c  gives  C  =  v/v  =  0.6.  Clearly,  larger  distances  A  should  result  in  smaller  effective  internal 
stresses  while  smaller  applied  stresses  should  give  increased  ratios  Oi/er.  The  overall  net  effect 
will  be  rather  complex  to  assess.  Here,  it  is  considered  that  under  the  conditions  described  the 
net  effect  on  the  overall  stress  exponent  of  the  creep  rate  will  be  to  elevate  it  from  n  =  3  for  the 
model  presented  above  to  n  =  (m  +  2)  =  (2.5  +  2)  =  4.5  to  n  =  (3  +  2)  =  5.  This  consideration 
and  the  factor  C,  when  incorporated  into  the  creep  model  of  smooth  climb  given  by  Eqn  (7) 
will  change  it  finally  to 


-  cos24>)3(2  -  3cos20) 
ft2  n  kT  ln{rs/rc) 


(10) 
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It  is  this  expression  that  is  compared  with  the  experimental  results  in  Section  5.4  below. 

5.4  Evaluation  of  the  creep  model 

The  final  creep  model  of  Eqn(10)  is  evaluated  now  and  compared  with  experimental  results.  In 
the  comparison  the  expression  is  evaluated  for  both  1400°C  and  1520°C.  The  following  model 
parameters  and  material  properties  are  used  in  the  comparison: 

H  =  =  73GPa  at  1400°C  and  63.4GPa  at  1520°C  [17] 

b  =  5.12  x  10-lom 
Q  =  1.06  x  10~29m3 
q  —  0.721  as  determined  in  Section  4.1 
(3  =  Q/b 3  =  0.079 
A  =  |(27r)3(g3//0)  =  7.85  x  102 
<j)  =  57.7° 

rs  =  0.5  —  2.0/im,  as  estimated  from  micrographs  (Fig.  7) 
rc  «  b 

D  =  10~14cm2/sec  at  1400°C  and  9  x  10_14cm2/sec  at  1520°C  determined  from  Oishi  and 
Kingery  directly  for  their  polycrystalline  material  [25]  for  which,  parenthetically  the  activation 
energy  was  Q=110kcal/mol,  considerably  larger  than  our  value. 

Finally,  we  take  n  =  4.5  as  suggested  from  the  analysis  in  Section  5.3  above,  and 
C  =  0.6. 

The  calculated  steady  state  creep  rate  relations  for  1400°C  and  1520°C  are  presented  in  Fig. 
6  as  the  broken  lines  in  comparison  with  the  experimental  data. 

Considering  the  several  uncertainties  in  the  model  and  material  parameters  the  agreement 
between  model  and  experimental  results  is  very  good. 

6  Discussion 

The  experimental  measurements  of  steady  state  creep  rates  in  the  Al203/c-Zr02  eutectic,  while 
few,  have  given  a  good  measure  of  the  response  of  this  material  in  the  temperature  range  of 
interest.  Taking  account  of  the  possible  interactions  of  the  c-Zr02  phase  with  the  Al203  we 
have  developed  a  steady  state  creep  model  based  on  climb  of  pyramidal  edge  dislocations  in 
Al203  alone.  In  the  model  it  is  proposed  that  the  c-Zr02  phase  of  relatively  large  size,  in 
the  micron  range  will  most  likely  undergo  nearly  complete  stress  relaxation  in  the  1400°C 
range,  resulting  in  an  initial  transient  creep  with  relatively  short  time  constant  of  roughly 
200  seconds  [2].  On  the  other  hand  the  sub-micron  size  domains  of  the  c-Zr02  will  most 
likely  be  too  small  to  undergo  independent  deformation  processes  by  crystal  plasticity  and  will 
remain  non-deformable  to  impede  the  climb  motion  of  the  pyramidal  system  dislocations  in  the 
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AI2O3  much  like  dispersoids.  The  observed  creep  rate  stress  exponents  in  the  range  of  4-5  are 
indicative  of  this  behavior.  In  this  sense  it  would  appear  that  the  eutectic  through  its  unique 
morphology  should  be  more  creep  resistant  than  sapphire  of  [0001]  axis  orientation.  To  explore 
this  comparison  the  creep  model  is  re-interpreted  for  application  to  sapphire  in  which  the  climb 
motion  of  the  dislocations  is  expected  to  be  quasi-smooth.  Such  a  modification  would  give  a 
steady  state  creep  expression  of 


•  _  ~  2  -  3  cos2<ft) 

tsaP  ln(rs/rc) 


(11) 


where  A'  =  §(27 r)3//?  =  2.09  x  103  (for  q=1.0) 

and  D1,  the  diffusion  constant,  must  be  chosen  for  sapphire  single  crystals  not  having  the 
benefit  of  diffusion  short  circuits  along  interfaces.  Prom  Oishi  and  Kingery  [25]  we  obtain 
D'  =  3.5  x  10-17cm2/sec  at  1400°C,  or  a  factor  of  3.5  x  10-3  lower  than  for  polycrystalline 
material  or  for  the  eutectic  with  a  large  volume  concentration  of  interfaces.  For  these  alterations, 
but  for  all  other  factors  remaining  the  same  a  steady  state  creep  expression  can  be  determined 
for  sapphire  crystals  of  [0001]  orientation  at  1400°C  which  is  shown  as  the  dashed  line  in  Fig. 
6,  confirming  the  expectation  of  a  lower  creep  resistance  in  comparison  with  the  Al203/c-Zr02 
eutectic. 


The  creep  model  that  was  presented  is  based  on  the  existing  evidence  of  the  measured  creep 
rates,  their  stress  and  temperature  dependence  and  on  the  findings  of  the  earlier  work  of  Fire¬ 
stone  and  Heuer  [3]  on  single  crystal  sapphire.  However,  the  work  leaves  many  unanswered 
questions.  Foremost  among  these  are  the  details  of  the  continuance  of  the  all-important  fluxes 
of  climbing  dislocations,  their  origins,  their  form  of  maintenance  and  the  fine  structure  of  cores 
of  dislocations  such  as  details  of  jogs  where  the  actual  climb  steps  occur  by  O  ion  vacancy  emis¬ 
sion.  Preliminary  energetic  considerations  of  dislocation  emission  from  ledges  on  the  ubiquitous 
interfaces  or  from  misfit  dislocations  along  interfaces  have  indicated  very  large  energy  barriers 
under  the  prevailing  local  stresses.  However,  other  forms  of  heterogeneous  nucleation  of  disloca¬ 
tions  from  interfaces  must  still  be  considered.  Other  possibilities  for  inelastic  behavior  such  as 
Coble  creep  were  discounted  because  of  the  apparent  general  absence  of  transverse  boundaries 
in  the  topologically  continuous  A1203  component  and  the  generally  linear  stress  dependence  of 
creep  flow  by  this  mechanism,  which  if  present,  appears  to  be  swamped  by  dislocation  climb 
flow.  Some  of  these  possible  mechanisms  as  well  as  the  missing  details  of  the  maintenance  of 
climb  fluxes  could  come  from  TEM  studies  which  were  in  progress  at  the  preparation  of  this 
final  report. 


7  Conclusions 

1.  In  the  eutectic  of  Al203/c-Zr02(Y203)  the  majority  phase  of  A1203  is  topologically  con¬ 
tinuous,  has  a  nearly  perfect  texture  of  [0001]  parallel  to  the  growth  direction  and  en- 
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capsulates  everywhere  the  minority  c-ZrC>2  phase  which  itself  has  a  growth  texture  of 
<  112  >. 

2.  The  interfaces  separating  the  phases  are  well  structured  but  are  incoherent  over  most  of 
the  area. 

3.  A  large  fraction  of  the  c-ZrC>2  phase  has  a  sub-micron  oriented  fibrillar  or  plate-like 
morphology  which  often,  but  not  always,  aggregates  into  colonies.  The  remaining  fraction, 
usually  surrounding  the  colonies  of  oriented  fibrils  or  platelets  is  of  a  coarser,  micron  size. 

4.  Upon  first  application  of  stress,  transient  creep  is  observed  which  is  attributed  to  stress 
relaxation  in  the  coarser  fraction  of  the  c-ZrC>2.  The  sub-micron  fraction  is  considered  to 
be  too  small  to  undergo  deformation  by  crystal  plasticity  and  is  expected  to  be  dormant. 

5.  In  steady  state  creep,  stress  changes  do  not  produce  additional  transients  within  the 
resolution  of  stress  measurements. 

6.  The  stress  dependence  of  steady  state  creep  rate  is  of  power-law  form  with  an  exponent 
in  the  range  between  4.5  to  5.0.  The  activation  energy  of  the  creep  rate  is  71.1  kcal/mol 
in  the  1200-1520°C  range  and  is  attributed  to  oxygen  ion  diffusion  through  the  AI2O3  and 
along  the  interfaces. 

7.  Because  of  the  [0001]  texture  of  the  AI2O3  phase  which  must  control  the  overall  creep 
rate  for  topological  reasons,  the  proposed  steady  state  model  is  based  on  the  climb  of  the 
(1/3)  <  1101  >  pyramidal  dislocations  in  AI2O3  since  the  latter  are  sessile  in  glide. 

8.  In  the  creep  model  which  is  a  generallization  of  the  steady  state  diffusional  creep  model 
of  Nabarro  [19],  the  stress  exponents  of  larger  than  3  are  attributed  to  the  necessity  of  re¬ 
peated  bowing  and  straightening  of  the  climbing  dislocations  in  AI2O3  as  they  circumvent 
the  sub-micron  c-ZrC>2  domains. 

9.  The  predictions  of  the  creep  model,  which  has  no  adjustable  constants,  agree  quite  well 
with  experimental  observations. 

10.  Finally,  a  mechanistic  comparison  of  the  creep  resistance  of  the  Al203/c-Zr02(Y203) 
eutectic  with  that  of  [0001]  sapphire  single  crystal  indicates  that  the  former  appears 
superior  because  of  the  obstructions  of  the  dispersoid-like  c-ZrC>2  domains  to  the  the 
climbing  dislocations  of  the  Al203  phase,  even  though  the  diffusion  constant  of  sapphire 
single  crystals  is  three  orders  of  magnitude  smaller  than  that  of  the  eutectic. 
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Table  1:  Recorded  Steady  State  Strain  Rates 


Stress  (MPa) 

Strain  Rate  e(sample  #)[e  average  rate]  (sec  l) 

1200°C 

200 

1.97  x  lO"10  (142-18)  i 

1400°C 

200 

3.6  x  lO"9  (187-29) 

225 

2.26  x  10~9  (187-27-2);  8.16  x  10~9  (187-27-1)  [5.21  x  10"9] 

250 

5.53  x  10~9  (187-29) 

275 

9.95  x  10~9  (187-27-2);  1.10  x  10~8  (187-27-1)  [1.048  x  10~B] 

1520°C 

150 

1.57  x  10"9  (142-16);  4.40  x  10"9  (142-16)  [2.98  x  10"9] 

200 

1.20  x  10“8  (142-16);  1.60  x  10“8  (142-16)  [1.40  x  10“8] 
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Table  2:  Material  Properties  for  Linear  Elastic  FEM  Analysis 


AI2O3 

c-Zr(>2 

Young’s  Modulus 

190GPa  at  1400°C 

0.05MPa  at  1400°C  a 

Poisson’s  ratio 

0.3 

0.4999999  a 

Yield  Stress 

5.0GPa  b 

5.0GPa 

Coefficient  of  Thermal 
Expansion,  °C~1 

9.1xl0_e  (a)  9.9x10-®  (c) [28] 

12.9-10-6  [28] 

“These  valurs  were  chosen  in  the  FEM  code  result  in  a  negligible  shear  modulus  of  83GPa 
for  c-Zr02  to  obtain  a  state  of  near  complete  shear  relaxation  while  leaving  the  bulk  behavior 
intact. 

6These  high  yield  stresses  were  chosen  in  the  FEM  code  to  guarantee  pure  elastic  response. 
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Figure  1:  An  SEM  micrograph  of  a  typical  morphology  of  the  Al203/c-Zr02  eutectic,  as  revealed 
in  a  transverse  section.  The  black  component  is  AI2O3.  Much  of  the  c-Zr02  component  (the 
bright  regions)  is  of  low  sub-micron  size  and  is  expected  not  to  deform  by  any  form  of  crystal 
plasticity. 
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Figure  2:  Sketch  showing  the  textural  arrangement  of  the  AI2O3  component  with  a  preferred 
[0001]  axis  and  c-ZrC>2  with  a  preferred  <  112  >  axis  found  by  electron  diffraction  analysis. 
Other  textures  for  c-Zr02  have  also  been  reported. 
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Figure  3:  A  cut-out  view  of  the  hot  zone  of  the  high  temperature  chamber  for  creep  experiments 
showing  also  the  cold  grips  of  specimens,  and  the  water  cooled  Cu  baffles  between  the  hot  zone 
and  the  grips. 


xIO'3 


Figure  4:  A  typical  strain  time  plot  for  steady  state  creep  at  two  different  stress  levels  in  an 
Al203/c-ZrC>2  eutectic  rod  at  1520°C 
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Figure  5:  A  typical  creep  strain  transient  in  an  A^Os/c-ZrC^  eutectic  under  a  stress  of  200 
MPa  at  a  temperature  of  1200°C.  The  horizontal  line  gives  an  estimate  of  creep  strain  increment 
due  to  stress  relaxation  in  the  coarse  fraction  of  c-ZrC>2. 
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Figure  6:  Plots  of  steady  state  creep  rates  as  a  function  of  stress  for  temperatures  of  1200, 
1400  and  1520°C  compared  with  measurements  reported  by  Sayir  and  Farmer  [2].  Broken  lines 
represent  predictions  of  creep  model  for  Al203/c-Zr02  at  1400°C  and  1520°C,  together  with 
model  prediction  for  sapphire  at  1400°C 
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Figure  7:  Determination  of  activation  energy  of  steady  state  creep  from  temperature  dependence 
of  steady  state  creep  rate. 


(a) 


(b) 


Figure  8:  A  Vickers  indentation  (a),  and  one  of  the  cracks  extending  from  one  corner  of  the 
indentation  (b).  The  crack  contour  is  used  to  estimate  the  fracture  toughness  of  the  eutectic. 
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(b) 


Figure  9:  An  internal  pore  (a),  and  a  surface  flaw  (b),  observed  on  the  premature  fractured 
surface  of  a  sample.  The  size  and  shape  of  these  defects  are  used  to  estimate  the  fracture 
toughness  of  the  eutectic 
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Figure  10:  An  AFM  fracture  surface  profile  characterization  (a),  and  corresponding  SEM  pic¬ 
ture  (b). 
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Figure  11:  Sketches  of  an  idealization  of  the  c-ZrC>2  domains  as  ellipsoidal  cylinders  in  the 
topologically  continuous  AI2O3  component,  for  the  purpose  of  a  FEM  study  of  the  internal 
stress  distribution  in  the  AI2O3  when  a  tensile  stress  az  is  applied. 
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Figure  12:  Stress  distribution  results  of  the  FEM  study;  a)  deviatoric  (Mises)  stresses  in 
two  components  of  the  eutectic;  b)  the  distribution  of  the  climb  stress  <7\  in  the  AI2O3. 


S,  Mises 

(Ave.  CriC.s  75%) 


-  +7 . 792e+02 
t +7 . 143e+02 

-  +6 . 493e+02 

-  +5 . 844  e+02 

-  +5 . 1956+02 

-  +4 . 545e+02 
. B96e+02 

-  +3 . 247e+02 
- +2 . 597e+02 

-  +1 . 946e+02 

-  +1.299e+02 
- +6.493e+01 
-+2.501S-04 


(a) 


Mean  Normal  Stress 
(Ave.  Crit.t  75%) 
+1.920e+02 
+1.4436+02 
+9.6676+01 
+4 . 900e+01 
+1 . 333e+00 
-4 . 633e+01 
-9.400e+01 
-1.417e+02 
-1 . 693e+02 
-2 . 370e+02 
-2 . B47e+02 
-3 . 323e+02 
-3 . 800e+02 


(b) 


Figure  13:  The  deviatoric  von  Mises  stresses,  (a);  and  the  mean  normal  stresses,  (b)  of  a 
temperature  reduction  from  1875°C  to  1400°C 
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Figure  14:  The  cation  sub-lattice  of  sapphire  showing  the  three  prominent  slip  system  in  the 
basal  plane;  in  the  prism  planes(broken  lines);  and  the  pyramidal  planes  (dashed)  lines. 
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Figure  15:  The  idealized  arrangement  of  ions  viewing  the  base  plane  from  above  [0001].  The 
large  open  circles  represent  an  underlying  layer  of  oxygen  anions;  the  empty  interstial  sites 
are  occupied  by  the  oxygen  anion  layer  immediately  above;  the  small  filled  circles  represent 
aluminum  cations  (in  the  real  crystal  these  do  not  lie  in  a  single  plane  but  are  alternately  asso¬ 
ciated  closer  to  one  or  the  other  nearest  oxygen  layers  where  they  produce  in-plane  distortions); 
the  empty  small  circles  represent  aluminum  holes  which  specify  the  scale  of  prefect  dislocation 
Burgers  vectors.  The  reference  rectangle  ABCD  represents  the  base  plane  of  the  simulation  cell 
for  the  base  plane  edge  dislocation,  while  the  rectangle  AA’CC’  represents  a  symmetry  plane 
in  the  simulation  cell  for  the  pyramidal  edge  dislocation  (after  Kronberg  1957). 
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Figure  16:  Pyramidal  plane  edge  dislocations  viewed  in  the  AI2O3  structure  from  the  [1120] 
direction.  Dislocations  are  on  (1102)  planes  with  b  =  1/3  <  1 101  >,  a)Unrelaxed  dislocation 
viewed  in  the  cation  (Al)  lattice  has  high  energy  but  could  glide  on  cation  ’’hole”  plane;  b) 
relaxed  dislocation  of  lower  energy  splits  core  into  two  fractional  edge  dislocations  lying  on 
parallel  cation  planes  of  very  high  resistance  making  dislocation  immobile. 


Figure  17:  Sketch  of  a  round  bar  of  a  A^Oa/c-ZrC^  eutectic  with  z  axis  parallel  to  the  [0001] 
direction  of  the  AI2O3.  The  plane  at  angle  <p  outlines  a  pyramidal  glide  plane.  A  local  tensile 
stress  o'  1  will  make  a  positive  edge  dislocation  climb  in  the  negative  X2  direction. 
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Climbing  Dislocation 


(C) 

Figure  18:  Sketch  depicting  the  required  contortions  that  a  climbing  edge  dislocation  needs 
to  suffer  in  bowing  through  the  gaps  between  the  c-Zr02  domains  that  we  consider  to  be 
impenetrable:  a)  critical  climb  configuration  to  bow  through  the  gap  between  two  c-Zr02 
dispersoids  at  a  spacing  A;  b)  at  the  point  when  the  climbing  dislocation  is  just  about  to  pinch 
off;  c)  the  cusped  dislocation  line  straightening  out  under  the  application  of  line  tension. 
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(b)  (c) 


Figure  19:  Consequence  of  the  repeated  bowing  out  and  straightening  of  the  climbing  dislocation 
considered  as  a  set  of  internal  resistance  a*  alternately  retarding  and  speeding  up  the  climbing 
dislocation:  a)  retardation  and  speeding  up  considered  as  unfavorable  and  favorable  internal 
stresses  cq  in  the  presence  of  a  climb  stress  cr;  b)  effect  of  Oija  on  the  overall  stress  exponent 
m  of  the  dislocation  velocity;  c)  effect  of  Oi/a  on  the  attenuation  factor  C  on  the  average  climb 
velocity,  (figures  reproduced  from  Li  [20],  courtesy  of  J.  Wiley  &  Sons) 
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